The dependence of the heUcity bound of force-free magnetic fields 

on boundary conditions 



This paper follows up on a previous study showing that in an open atmo- 
sphere such as the solar corona the total magnetic helicity of a force-free field 
must be bounded and the accumulation of magnetic helicity in excess of its upper 
bound would initiate a non-equilibrium situation resulting in an expulsion such 
as a coronal mass ejection (CME). In the current paper, we investigate the de- 
pendence of the helicity bound on the boundary condition for several families of 
nonlinear force-free fields. Our calculation shows that the magnitude of the he- 
licity upper bound of force-free fields is non-trivially dependent on the boundary 
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^ I condition. Fields with a multipolar boundary condition can have a helicity up- 

per bound ten times smaller than those with a dipolar boundary condition when 
^ helicity values are normalized by the square of their respective surface poloidal 

' fluxes. This suggests that a coronal magnetic field may erupt into a CME when 

^ I the applicable helicity bound falls below the already accumulated helicity as the 

result of a slowly changing boundary condition. Our calculation also shows that a 
' monotonic accumulation of magnetic helicity can lead to the formation of a mag- 

^ ■ netic flux rope applicable to kink instability. This suggests that CME initiations 

O ■ by exceeding helicity bound and by kink instability can both be the consequences 

of helicity accumulation in the corona. Our study gives insights into the observed 
^ ■ associations of CMEs with the magnetic features at their solar surface origins. 
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1. Introduction 

Magnetic helicity is a physical quantity that measures the topological complexity of a 
magnetic field, such as the degree of linkage and/or twistedness in the field (Moffatt 1985, 
Berger & Field 1984). In a previous paper (Zhang et al. 2006, hereafter referred to as ZFL) 
we proposed that in an open atmosphere such as the solar corona there is an upper bound on 
the total magnetic helicity that a force-free field can contain. The accumulation of magnetic 
helicity in excess of this upper bound would initiate a non-equilibrium situation, resulting 
in a coronal mass ejection (CME) as a natural product of coronal evolution. 

Our approach (Zhang & Low 2003, 2005, ZFL) shifts the traditional focus on the mech- 
anism for CME evolution from the storage of magnetic energy to the accumulation of mag- 
netic helicity, although the two types of considerations are not necessarily exclusive with 
each other. The advantage of using magnetic helicity as a more fundamental physical quan- 
tity is that from observations we know that the magnetic fields are emerging from the solar 
interior with a preferred helicity sign in each solar hemisphere (Pevtsov et al. 1995, Rust 
& Kumar 1996, Bao & Zhang 1998, Zhang 2006). As a result, an accumulation of the total 
magnetic helicity in the corona becomes unavoidable because the total magnetic helicity 
is approximately conserved in the corona during coronal processes including fast magnetic 
reconnection (Berger 1984). 

In this paper, we study how the magnitude of the upper bound of the total magnetic 
helicity depends on the boundary condition. Section 2 presents the model with two new 
boundary conditions. Results and analysis are presented in Section 3. In Section 4, a brief 
summary of the paper is given. 



The Model 



2.1. The governing equation 



Following Flyer et al. (2004) and ZFL, we use the famihes of power-law axisymmetric 
force-free fields to understand the basic physical properties of interest. 



With axisymmetry, the solenoidal magnetic field B in r > 1 can be written in the form 



of 



B = 



r sin 9 



idA _ dA ■ 

r 86 ^ (9r ' 



(1) 



where the fiux function A defines the poloidal magnetic field and the function Q defines the 
toroidal (or azimuthal) field. 
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Q is defined as a strict power-law in A with the form 

where n is an odd constant index required to be no less than 5 in order for the field to possess 
finite magnetic energy in r > 1 and 7 is a free parameter which we choose to be positive 
without loss of generality. This form of Q reduces the force-free condition to the following 
governing equation for the flux function A: 

^ + ^^^^ + 7A"-0 (3) 

This governing equation was solved numerically as a boundary value problem within 
domain r > 1 in Flyer et al. (2004) and ZFL, subject to the prescribed boundary flux 
distribution of 

A\r=i = sin^e . (4) 

This boundary flux distribution and its associated normal field distribution are plotted in 
the top panels of Figure 1. All force- free fields discussed in Flyer et al. (2004) and ZFL 
share the boundary flux distribution given by (jlj), which is the same as that for a dipole 
potential field corresponding to the solution of equation ([3]) with 7 = 0. The solutions of ([3]) 
with boundary condition (jll) shall be referred to as dipolar force-free fields or dipolar fields 
for short. We refer interested readers to Flyer et al. (2004) and ZFL for various properties 
of dipolar force-free fields for the cases of n = 5, 7, 9. 



2.2. Two new boundary conditions 

In this paper, we solve the same governing equation ([3]) using the numerical methods 
described in Flyer et al. 2004, 2005 but subject to two new and distinctively different 
boundary conditions. We use these numerical solutions to investigate how properties we 
discussed in ZFL would change with the new boundary conditions. 

The first new boundary condition has the flux distribution of 

A\r=i = sin^H . (5) 

This flux distribution and its associated normal field distribution are plotted in the middle 
panels of Figure 1. We see that this new boundary condition has its flux concentrated nearer 
to the equator than that of the dipolar field. This makes it more like solar active regions 
with its normal field strength much higher at equatorial regions than that near the poles. 
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We shall refer this family of power-law axisymmetric force-free fields as bipolar force-free 
fields or bipolar fields hereafter. 

The second new family of solutions are those that we shall refer to as multipolar force- 
free fields or multipolar fields. They are also the solutions to ([3]) but subject to the following 
boundary condition: 

A|^=i = sin^e X (1 -5 X cos^^) . (6) 

Its flux distribution and associated normal field distribution are plotted in the bottom panels 
of Figure 1. We see that fields with this boundary condition have both positive and negative 
magnetic fields in one hemisphere, distinctively different from those in dipolar force-free 
fields and bipolar force-free fields where fields in one hemisphere have the same magnetic 
polarity. 

Table 1 presents a comparison of the three boundary conditions. We see that their com- 
mon feature is that they all have A\r=i = at solar poles and A\r=i = 1 at the equator. The 
dipolar boundary and bipolar boundary differ by the contrast of their respective maximum 
Br\r=i values to their Br\r=i values at the northern pole. The multipolar boundary condition 
differs from the dipolar and bipolar boundary conditions by the existence of both magnetic 
polarities in each hemisphere, indicated by its negative minimum A value. 

In Table 1 we also present the total surface poloidal flux (Fp) for the three boundary 
conditions, where 

TT 

Fp= I Br{>0)ds = 271 \Br\smedd . (7) 

Jr=l Jo 

We see that dipolar and bipolar fields have the same total surface poloidal flux {2%) and the 
multipolar fields have a larger total surface poloidal flux {5.2tt). In the latter development 
of the paper, we will normalize calculated magnetic helicity by in order to make fields 
with different boundary conditions comparable. 

3. Results and Analysis 

As in ZFL, for each new boundary condition, we numerically solve ([2]) for three cases: 
n = 5, n = 7 and n = 9. In each case, the numerical method, that is, the Newton's iteration 
combined with a pseudo-arc length continuation scheme, guarantees the completeness of each 
solution branch generated by the 7 values. 

Also as in ZFL, for each obtained solution, we calculate three physical quantities of the 
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field. They are: total magnetic energy 



total azimuthal flux 



ai au = \ I — 

y m 

and total relative magnetic lielicity 



V>i y m + 1 J^yi sm 9 



7 f ^^m+2drdd 

ryl V + 1 Jr>l sin 9 



Hr = Ati AB^rdrde = AnJ — — / . (10) 



The derivation of these formula can be found in ZFL. The only difference is that since we 
are also considering multipolar fields the absolute value of B^. (that is, \Br\) is introduced 
when calculating total azimuthal flux. 

In the geometric simplicity of these force-free fields, the equilibrium in each case is due 
to the magnetic tension force of the poloidal flux confining the magnetic pressure of the 
azimuthal flux. The magnitude of the tension force is sensitively related to the poloidal flux 
at the inner boundary (r = 1) that serves as an anchor agent. Moreover, the poloidal flux 
and its tension force become weaker in the outward radial direction. Hence, if the azimuthal 
flux becomes too large it can not be confined by such a mechanism as was shown in ZFL 
for dipolar fields. The magnetic pressure is independent of the sign of the field. Therefore, 
the total unsigned azimuthal flux given by is more relevant for the consideration of flux 
confinement than the total signed azimuthal flux. 



3.1. Helicity upper bound of bipolar fields 

Figure 2 presents the variations of the total magnetic energy {E), total azimuthal flux 
(F^) and total magnetic helicity (Hr) versus 7 along the solution curve for bipolar force-free 
flelds with n = 5, n = 7 and n = 9. The flgure is similar to the Figure 2 of ZFL except that 
Figure 2 in ZFL is for dipolar force-free flelds. Each point along the solution curve, denoted 
by a plus symbol in the flgure, represents a solution to (j3]). By solution curve we mean that 
all solutions along the curve are obtained with the same boundary condition and the same 
constant index n but with a monotonically increasing magnitude of total azimuthal flux. 

From these curves of solutions we see that there may be upper bounds on the total 
magnetic energy, total azimuthal flux as well as total magnetic helicity for bipolar flelds, as 
that we have suggested in ZFL for dipolar flelds. 
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Prom Figure 2 we can also see that the magnitude of the hehcity upper bound for bipolar 
fields is smaller than that for dipolar fields. While the upper bound for dipolar fields is close 
to 14, the upper bound for bipolar fields only approaches 9. To illustrate this further, we 
plot in Figure 3 the variation of the total magnetic helicity {Hji) versus azimuthal flux (F<^) 
along the solution curve for fields with the dipolar (top panels) and bipolar (bottom panels) 
boundary conditions. Here we have normalized the values of total magnetic helicity (Hji) 
of each field by the square of their corresponding surface poloidal fluxes (F^ ) as those in 
Demouhn et al. (2002), van Driel-Gesztelyi et al. (2003) and Demoulin (2007). We see 
that while the upper bound of Hr/F^ is about 0.35 for dipolar fields, for bipolar fields it is 
significantly lower at 0.22. 

This result suggests a dependence of the helicity upper bound on the boundary condition, 
which in our view gives an insight into observed associations between CMEs and magnetic 
features at their solar surface origins. Observationally we know that CMEs can be triggered 
by flux emergence (e.g. Feynman & Martin 1995, Subramanian & Dere 2001) and converging 
motions (e.g. Martin 1990). Different theoretical models have also proposed that CME-type 
eruptions can be triggered by various surface field variations (e.g. Chen & Shibata 2000, 
Amari et al. 2000, 2003a, 2003b). 

We suggest that the common physics underlying the different mechanisms associated 
with such surface field variations is the dependence of helicity upper bound on the boundary 
condition. When a magnetic field has accumulated a certain amount of magnetic helicity 
(but not yet enough for an eruption) then a change of the boundary condition could lower 
the helicity upper bound, resulting in a non-equilibrium situation and hence a CME eruption 
under the new boundary condition. For example, if Hr/ were 0.3 for a dipolar boundary 
condition, then an evolutionary change to a bipolar boundary condition would result in a 
CME eruption because the applicable upper bound on the conserved total hehcity has been 
reduced as suggested by our numerical experiments. 

A note to address here is that although in this paper we have emphasized the role of 
boundary condition variations this does not mean that the role of magnetic helicity accumu- 
lation becomes less important. A change of the boundary condition may bring in an eruption 
only when the field has accumulated enough helicity for an eruption under the new boTind- 
ary condition. If not, the field does not erupt even when the boundary fiux distribution is 
changing. This is consistent with the observation (Zhang et al. 2007) that although fiux 
emergences are indeed found to be associated with CME eruptions, the same rate of fiux 
emergence can also be found when there is no CME or solar activities. This means that fiux 
emergence may be a trigger of a CME eruption, but flux emergence alone do not guarantee 
an eruption. 
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Another interesting result from our calculations is that these normahzed Hu/Fp helici- 
ties, estimated from simple axisymmetric power-law force-free fields, lie close to those esti- 
mated from observations. Van Driel-Gesztelyi et al. (2003) and Dcmoulin (2007) pointed out 
that the H^/ numbers, estimated from the extrapolated magnetic fields based on observed 
photospheric magnetograms, are between 0.02 to 0.2. Our helicity upper bound numbers of 
dipolar and bipolar fields are just a little higher. Notice that the numbers estimated from 
the observations may be somehow underestimated because of the limited spatial resolution 
of the observed magnetograms. So there seems to be consistency between the theoretical 
Hn/Fp hehcities and those estimated from observations. 

Figure 4 presents four field configurations selected from the n = 9 solution curve, po- 
sitions of which along the curve are illustrated in Figure 2. We see that starting from the 
potential field (Panel A) the curve first reaches the maximum-7 field (Panel B) and then 
the field with maximum total magnetic energy (Panel C) and then the one with maximum 
total azimuthal fiux (Panel D). A clear bubble (representing a fiux rope) is presented in the 
field of Panel C but not in the field of Panel D where the latter actually possesses more 
total azimuthal fiux than the former. This tells us that although the existence of a fiux rope 
(or fiux ropes) in the low corona does represent a storage of a certain amount of magnetic 
helicity (see more discussions in Zhang & Low 2003) it is not necessary that they are present 
in the field with a maximum helicity storage. 



3.2. Helicity upper bound of multipolar fields 

As in Figure 2, we present in Figure 5 the variations of the total magnetic energy, total 
azimuthal flux and total magnetic helicity versus 7 along the solution curve with n = 5, 
n = 7 and n = 9 but for the multipolar fields. We see that these curves also suggest the 
existence of upper bounds on the total magnetic energy, total azimuthal flux as well as total 
magnetic hehcity as those for dipolar flelds and bipolar flelds. 

The flgurc also shows that the magnitude of the helicity upper bound for multipolar 
flelds is smaller than that for dipolar flelds. This reduction is even more evident in Figure 
6 where we plot Hr/ F^ versus F<^. We see that the upper bound of Hr/ F^ for multipolar 
flelds is below 0.04, almost ten times smaller than that of dipolar flelds. 

Such a severe reduction of Hr/F^ upper bound not only further conflrms our previous 
result that the helicity upper bound is dependent on the boundary condition, but also brings 
our theoretical Hr/ F^ value even closer to those estimated from observations (Regnier et 
al. 2005). Furthermore, the severe reduction of the helicity upper bound in terms of Hji/Fp 
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values may also explain why solar eruptions such as CMEs are more likely to happen in 
complicated active regions where the multipolar field, by the above property, will take less 
time to reach its helicity bound, producing an eruption. 

As in Figure 4, Figure 7 presents four field configurations selected from the n — 7 
solution curve of the multipolar fields. We see again a clear bubble in the field of Panel 
C but not in Panel D. Therefore, as with bipolar fields, multipolar fields with maximum 
helicity storage need not contain a flux rope. 

3.3. Kink instability 

A rope of highly helical field is susceptible to an instabihty that causes the rope to kink 
(Friedberg 1987). From elementary calculations, this kink instability sets in if a critical twist 
is exceeded (T > T,,). The exact value of T^. depends on the detailed field models, and could 
increase from the traditional Kruskal-Shafranov limit — 27r to — 2.57r (Hood & Priest 
1981) and % = 4.87r (Mikic et al. 1990). 

Since we have helical flux tubes (or flux ropes) present in our solutions, it is interesting 
to investigate whether these flux ropes have exceeded the kink instability. Figure 8 presents 
the variation of the average twist (T) versus for two flelds. One of the flelds is the bipolar 
n = 9 maximum-energy fleld, the one presented in Panel C of Figure 4. The other is the 
multipolar n = 7 maximum-energy fleld, the one presented in Panel C of Figure 7. Here 6o is 
the angle from the equator. The average twist (T) is obtained from H'j^/{Fp)'^, where H'j^ is 
the total relative magnetic helicity in domain Jl', enclosed by the r — 1 surface and the flux 
surface with A — A\j.=i^0=0q, and is the total surface {r — 1) polodial flux of the domain 

n'. 

We see that in both flelds the average twist of the central part of the fleld (that is, 
where the flux rope is located) has exceeded the kink instabihty criteria, Tc — 2.57r of Hood 
& Priest (1981) or = 4.87r of Mikic et al. (1990). This tells us that with the accumulation 
of a certain amount of magnetic helicity, the flux rope formed in the fleld can possess a 
twist number that is larger than the kink instability criteria. If other necessary conditions 
are favorable, for example, if the field has accumulated enough free magnetic energy, an 
eruption may happen even before the field has reached its helicity upper bound state. In 
that sense, reaching the helicity upper bound state may not be a necessary condition for 
eruption, but the helicity upper bound is a sufficient condition upon which an eruption will 
become unavoidable. This also shows that CME eruptions that are initiated by the kink 
instability (e.g. Torok and Kliem 2005, Fan & Gibson 2007) or by the existence of helicity 
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upper bound could both be viewed as the consequences of magnetic hehcity accumulation 
and they are not mutually exclusive. 



4. Conclusion 

In this paper, we continue our study on the hydromagnetic origin of CMEs in terms of 
magnetic helicity accumulation. As in a previous paper (ZFL), we numerically solve ([3]) to 
get families of axisymmetric power-law force-free fields, but subject to two new boundary 
conditions. 

By analyzing and comparing obtained solutions for three different boundary conditions 
we conclude the following: 

1. The suggestion that there may be an upper bound on the total magnetic helicity for 
force-free fields is also found for the two new boundary conditions. 

2. The magnitude of the helicity upper bound of force- free fields is non-trivially dependent 
on the boundary condition. In our examples, the fields with a surface flux distribution 
more like a simple active region (bipolar fields) have their helicity upper bound smaller 
than that of fields with dipolar boundary condition. For multipolar fields, the helicity 
upper bound {Hr/ F^) can be ten times smaller than that of dipolar fields. These results 
provide some insights into the observed association of CMEs with flux emergence and 
surface field variation. These results also suggest a physical reason why eruptions are 
more likely to happen in complicated active regions. 

3. CME initiations by kink instability and by the existence of a helicity upper bound 
can both be the result of magnetic helicity accumulation in the corona. They do not 
exclude each other. 
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Table 1: Comparison of the three boundary conditions 
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Fig. 1. — Flux distribution (left panels) and normal field distribution (right panels) of dipolar 
fields (top panels), bipolar fields (middle panels) and multipolar fields (bottom panels) at 
r = 1. 
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Fig. 2. — Variation of the total magnetic energy (E), azimuthal flux (F^) and total magnetic 
helicity (Hr) vs 7 along the solution curve for = 5 (top panels), n = 7 (middle panels) and 
n = 9 (bottom panels) fields with the bipolar boundary condition. Each point, denoted by a 
plus symbol in the plot, represents a solution to (I3l). Letters (A, B, C and D) in the bottom 
panels indicate the positions along the solution curve of the four fields plotted in Figure 4. 
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Fig. 3. — Variation of the total magnetic helicity (Hji) vs azimuthal flux (F^) along the 
solution curve for fields with the dipolar (top panels) and bipolar (bottom panels) boundary 
conditions. Here the total magnetic helicity (Hr) of each field has been normalized by the 
square of their corresponding poloidal flux (F^). 
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Fig. 4. — Field configurations of four fields selected from the n — 9 bipolar fields. These 
contours of flux function A represent the lines of force of the axisymmetric field projected on 
the r — 9 plane. Panel A: the potential field; Panel B: the field with the maximum 7 value 
among all n = 9 solutions; Panel C: the field possessing the maximum total magnetic energy 
among all n = 9 solutions. The thick line in this plot outlines the range within which the 
average twist is greater An (see text in Section 3.3 for details); Panel D: the field possessing 
the maximum azimuthal flux among all n = 9 solutions. 




Fig. 5. — Same as Figure 2, but for fields with the multipolar boundary condition. 
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Fig. 6. — Variation of the total magnetic helicity (H^), normalized by the square of poloidal 
flux (Fp), vs azimuthal flux {F^) for fields with the multipolar boundary condition. 
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Fig. 8. — Variation of average twist (T) vs 6*0 (angle from the equator) for the bipolar n — 9 
maximum-energy field (top panel) and multipolar n — 7 maximum-energy field (bottom 
panel). See text for the calculation of T. The dashed lines indicate the critical twist of kink 
instability where = 2.57r (Hood & Priest 1981) or Tc = 4.87r (Mikic et al. 1990). See also 
Panel C in Figures 4 and 7 for field configurations. 



